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1 Introduction 

In the recent years, there has been some interest for a class stochastic processes which are 
meant to serve as models for the evolution of an object that breaks down into smaller pieces, 
randomly and repeatedly as time passes. We refer to the monograph [6] and the survey [5] 
for a detailed account and references. Several crucial hypotheses have to be made in order 
to deal with models that can be analized by standard probabilistic techniques. Typically, one 
assumes that the process enjoys the branching property, in the sense that the dynamics of a 
given fragment do not depend on the others. A further important assumption which is made 
in [6], is that each fragment is characterized by a real number which can be viewed as its size. 
The latter requirement does not allow us to consider geometrical properties like the shape of a 
fragment, although such notions could be relevant for describing how an object breaks down. 



Figure 1: Example of a fragmentation of a square into squares and triangles 



In the simpler case when time is discrete, one can analyze a fragmentation chain using the 
framework of branching random walks or that of multiplicative cascades. In this setting, it is 
therefore natural to enrich the model by assigning to each fragment a type, which, for instance, 
may describe its shape, and let the evolution of each fragment depend on its initial type. The 
study of the latter can then be developed directly by translating the literature on multitype 
branching random walks or cascades (see e.g. [21 El [9]). 

However, we shall be interested here in the much more delicate case where time is continuous 
and each fragment may splits immediately, a situation which cannot be handled directly by 
discrete techniques based on branching processes. In the monotype setting, Kingman's theory 
of exchangeable random partitions provides the key for the construction and the study of 
fragmentation processes in continuous time; this was pointed out first by Pitman [13], see 
Chapter 3 in [6] for a complete account. In the first part of this work (Sections 2 and 3), we 
shall briefly explain how Kingman's theory can be extended to the multitype setting (for any 
finite family of types), and how this extension enables us to develop an adequate framework 
for multitype fragmentation processes. In short, the main result states that the dynamics of 
a homogeneous multitype fragmentation are characterized by a family of erosion coefficients 
and a family of dislocation measures. Each erosion coefficient describes the rate at which a 
fragment with a given type melts down as time passes, and each dislocation measure specifies the 
statistics of its sudden splits. Once the correct setting is found, statements are straightforward 
modifications of that in the monotype situation, and for the sake of avoiding what would be 
essentially a lengthy and boring duplication of existing material, our presentation will be rather 
sketchy and proofs will be omitted. Non-specialist readers may wish to consult first Chapters 
2 and 3 of [S] for getting the flavor of the arguments. 

The second part of this work (Section 4) is devoted to the study of the tagged fragment, i.e. 
the fragment which contains a point which has been tagged at random and independently of the 
fragmentation process, and its applications. It departs more significantly from the monotype 
situation, in the sense that the evolution of the tagged fragment is now given in terms of a 
Markov additive process (instead of a subordinator), which depends explicitly on the charac- 
teristics of the fragmentation. The central limit theorem for Markov additive processes then 
enables us to determine the asymptotic behavior of certain multitype fragmentation processes, 
extending an old result of Kolmogorov [13J in this area. We will also develop the natural connex- 
ion with multitype branching random walks from which we derive some sharp large-deviation 
estimates based on the work of Biggins and Rahimzadeh Sani [9] . 
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Throughout this text, we shall consider a finite family of types, say with cardinal k + 1 > 2, 
which can thus be identified with {0, 1, ... , k}. The type is special and will only be used in 
peculiar situations. As it has been mentioned above, it may be convenient to think of a type 
as a geometrical shape (see Figure 1 above for an example), but the type can also be used, 
for instance, to distinguish between active and inactive fragments in a frozen fragmentation 
(see [10] for a closely related notion in the setting of coalescents). Last but not least, it was 
observed recently by Haas et al. [UJ that homogeneous fragmentions bear close connexions with 
certain continuum random trees, a class of random fractal spaces which has been introduced 
by Aldous. It is likely that more generally, multitype fragmentations can be used to construct 
some multifractal continuum random trees, following the analysis developed in [llj. 

2 Kingman's theory for partitions with types 

The purpose of this section is to provide a brief presentation of an extension of Kingman's 
theory (see [12] or Section 2.3.2 in [6]) to partitions with types, which is a key step in the 
analysis of random fragmentations. 

2.1 Partitions with types 

We shall deal with two natural notions of partitions with types, which correspond to two 
different points of view. The first one focuses on the masses (and the types) of the components, 
whereas, roughly speaking, the second one corresponds to a discretization of the object which 
breaks down. 

We call any numerical sequence x = (x\, x 2 , ■ ■ ■) with 

oo 

Xi > x 2 > ■ ■ ■ > and x n < 1 

l 

a mass-partition, and write V m for the space of mass-partitions. A mass-partition with types 
is a pair x = (x, i) with x = (x\, x 2 , . . .) G V m and i = (ii,i 2 , ■ ■ •) a sequence in {0, 1, ... , k}, 
such that for every n G N = {1, 2, . . .} 

x n = 0^i n = 0, (1) 

and further 

i.e. the sequence (x±, ii), (x 2 , i 2 ), . . . is non-increasing in the lexicographic order. We shall write 
indifferently 

x= (x, i) = ((x 1 ,i 1 ),(x 2 ,i 2 ),...) 

by a slight abuse of notation. 

We should think of x n as the size of the rz-th largest component of some object with total 
mass 1 which has been split, and of i n as its type. A component with size means that it is 
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absent or empty, and thus has the special type 0. Note that a mass-partition x can be improper, 
in the sense that x n < 1- Then the mass-defect xo = 1 — x n is called the mass of dust, 
where the dust is viewed as a set of infinitesimal particles. It may be convenient to think that 
the special type is also assigned to these infinitesimal particles. 

We write V m for the space of mass-partitions with types and endow it with the following 
distance. Let (ei, . . . , e&) denote the canonical basis of the Euclidean space IR fc , and associate 
to any mass-partition with types x £ V m the probability measure on the axes of the unit cube 

oo 

y?x = x 5 + ) y x n 5 Xnein , 

n=l 

where Xq = 1 — * s the mass °f dust. Then we define the distance cf(x, x') for every 

x, x' £ V m as the Prohorov distance between the probability measures ip^ and if^>, which 
makes (V m , d) a compact space. We stress that the distance d is strictly weaker than other 
perhaps simpler distances on V m such as, for instance, 

max{|x n - x' n \\ {in=i , n} + (x n + x' n )l {in ^ n} : n £ N} . (2) 

Next we turn our attention to the second notion of partition. We call any subset of N = 
{1, 2, . . .} a block. A partition of a block B C N is a sequence n = (tti, 7T2, . . .) of pairwise 
disjoint blocks with U7r n = B, which is ranked according to the increasing order of the least 
elements, i.e. inf n m < inf ir n whenever m < n (with the usual convention that inf = oo). 
We write Vb for the space of partitions of B. 

Given a partition it = (tti, 7r 2 , . . .) £ Vb, we can assign to each block n n a type i n £ {0, . . . , k}, 
with the following convention which is related to §B) : 

i n = -v^ 7r n is either empty or a singleton. (3) 

We write n = (ir, i) = (^2^2), • • •) and call rr a partition with types of B. We denote 

by Vb the space of partitions with types of some block B. 

For every block BCN and every partition tc = (tti, . . .) of N, we define k\b, the restriction 
of 7r to B, as the partition of B whose blocks are given by 7c n n B, n £ N. If W — (tt, i) is now 
a partition with types, we assign types to the blocks of the restricted partition h\b as follows. 
The type of 7r n D B coincides with the type i n of the block n n if n n n B is neither empty nor a 
singleton, and is otherwise in order to agree with (J3]). We then write tt\b for the restriction 
to B of the partition with types 7f. 

For every pair (ff , 7f') of partitions with types, we define 

d(7f, 7f') = 1/ SUp{n £ N : Tt\[ n ) = 7f|[ n ]} , 

x It may be worthy to point out that our choice for the distance is well-adapted to the requirement (TI]). 
Typically, denote for n G N by x*-™'^ the mass-partition with types which consists in n identical fragments with 
mass 1/n and fixed type i € {1, . . . , k}. Then xS n 'V converges as n — > 00 to the degenerate partition of pure 
dust (and type 0). Such a natural convergence would fail if we had chosen a stronger distance on P m like ((2])- 
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where [n] — {1, . . . , n} and 1/ sup N = 0. Note that, as the type assigned to singletons is always 
0, the identity 7f|[i] = ftL-, holds in all cases and thus dijr, 7f') < 1. It is easily checked that 
d(n, it') defines a distance which makes V?q a compact set; see Lemma 2.6 in [6] on its page 96. 
Finally, we say that a block BCN possesses an asymptotic frequency if and only if 

\B\ = lim n" 1 Card( J Bn [n]) 

?i— >oo 

exists. If all the blocks of n E P N possess an asymptotic frequency, then we say that 7t has 
asymptotic frequencies, and we write \lt\^ = (x, i) for the sequence of the asymptotic frequencies 
and types of the blocks of 7t ranked in the non-increasing lexicographic order. Note from Fatou's 
lemma that \n n \ < 1 is a mass-partition and thus G V m . 

2.2 Exchangeability and paintbox construction 

A finite permutation is a bijection a : N — » N such that u{n) = n when n is sufficiently 
large. The group of finite permutations acts naturally on the space of partitions with 
types. Specifically, we write a -1 for the finite permutation obtained as the inverse a. Given 
an arbitrary partition with types of N, 7f = (ir, i), cr _1 maps each block ix n of it into a block 
o" _1 (7i"n) of some partition denoted by (t(tt). We decide to assign the type i n of the block 7r n to 
the block a~ l {ii ri ). This way, we obtain a partition with types denoted by <t(tt). 

A measure on Vn is called exchangeable if it is invariant under the action of finite per- 
mutations. Following Kingman [12J, we can associate to every mass-partition with types 
x = (x, i) G V m an exchangeable probability measure on P N by the paintbox construction. 
Specifically, introduce a pair of random variables (£, r) with values in Z + X {0, 1, . . . , k} whose 
distribution is specified by the following : 

oo 

P((£, r) = (n, i n )) = x n for every neN and P((f , r) = (0, 0)) = 1 - ^ x n . 

n=l 

Then consider a sequence (^i, ti), . . . of i.i.d. copies of (£, r) and define a random partition with 
types 7r = (7r, i) by declaring that two distinct integers £, m are in the same block of it if and 
only if £ m = & > 1, and then decide that the type of that block is r m = T£. Integers £ such 
that ^£ = form the class of singletons of n, and their type is of course 0. Similarly, if some 
block of 7r is empty, then its type is necessarily by our convention. The distribution of 7f will 
be denoted by and called the paintbox based on x. 

A slight variation of this paintbox construction can be illustrated as follows. Suppose for 
simplicity that the mass-partition with types x can be represented by splitting some geometric 
object, for instance a rectangle with unit area, into smaller components, for instance squares, 
rectangles and triangles. Each component has an area and a shape which we called a type. 
Imagine that we pick at random a sequence of i.i.d. uniform points U\, U2, ... in the initial 
object. A random partition with types is obtained by declaring that two distinct indices are 
in the same block of the partition if the corresponding random points belong to the same 
component of the object, and the type of this block is then the type of this component. See 
Figure 2 below. 
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Figure 2 : paintbox with types for a partition of a rectangle; 
types : square = \, rectangle = 2 and triangle = 3 

7f| [n] = (({1, 4, 10}, 2), ({2, 6}, 1), ({3, 8}, 2), ({5, 7}, 3), ({9, 11}, 3), (0, 0), . . .) 

By the law of large numbers, for every positive integer £, the block B = {m G N : £ m = £} 
has an asymptotic frequency 

\B\ = lim n" 1 Card( J B n [n]) = x e , 

n— >oo 

and clearly the type i^. One can arrange the sequence of the pairs (asymptotic frenquency, 
type) of the blocs of 7f in the non-increasing lexicographic order and then it coincides with x. 

Another important observation is that |7fi|, the asymptotic frequency and the type of the 
first block tti of a paintbox based on a mass-partition with type x, has the distribution of a 
size-biased sample of x, that is 

oo 

Px(|tti| = (x n ,i n )) = x n for every n eN and Px(|7h| = (0,0)) = x = 1 - ^x n . (4) 

l 

Plainly, if a is a finite permutation, then (£o-(i), T a m), (£, a (2), r o-(2)), • • • is again a sequence of 
i.i.d. copies of (£, r) and the corresponding partition with types is given by (t{tt). Thus px 
is an exchangeable probability measure on T>n, and more generally any mixture of paintboxes 
produces an exchangeable probability measure on V m . The converse to the latter assertion is 
a slight variation of the fundamental theorem of Kingman [12], see e.g. Theorem 2.1 in [6] on 
its page 100. 
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Theorem 1 Let p be an exchangeable probability measure on V®. Then p-almost every fc G Vn 
possesses asymptotic frequencies, and if g stands for the distribution of \ti\^ under p, then there 
is the following disintegration of the measure p : 

p(dfc) = / px(dvf)^(dx) , ttgPn- (5) 

Conversely, for every probability measure g on V m , ((SI) defines an exchangeable probability 
measure on V®. 



We next turn our attention to an extension of Kingman's theorem to certain sigma-finite 
measures on P N . In this direction, it is convenient to denote for every type i G {1, . . . , k} and 
every block B that is neither empty nor a singleton, by 1 b,% the partition with type of B given 
by ((£?, i), (0, 0), . . .). We also write lj = ((1, i), (0, 0), . . .) G V m for a related mass-partition 
with types. For every n G N, we denote by e^j) for the partition with types of N which has 
exactly two non-empty blocks, (N\{n},i) and ({n},0). The exchangeable measure on V® 

nGN 

will be referred to as the erosion measure with type i. 

The following extension of Theorem [1] to certain possibly infinite measures is the multitype 
version of Theorem 3.1 in |6J on its page 127. Recall the notation tt\b for the partition with 
types restricted to some block B, and that [2] = {1,2}. 

Theorem 2 Fix a type i G {1, . . . , k} and let pi be an exchangeable measure on V® such that 

Ati({lN,i}) = and p { (n eVn ■ %] ^ ![2],i) < oo . (6) 
Then the following holds: 

(i) pi-almost every partition tt G Pn possesses asymptotic frequencies. 

(ii) Let be the image measure of pi by the mapping ft — > \7t\K The restriction 

i/<(dx) = l^ijl/iil^dx) (7) 

of \pi\ l to V m \{h} fulfills 

and there is the disintegration 

l{|*|Vij^( d7r ) = / Px(dvf)z/j(dx) . 

(iii) There is a real number Cj > such that 

l { | ft H =1 .}//i(d7f) = c iei (d7r). (9) 



Vn 



[1 - xil {il= j})^(dx) < 00 
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Conversely, for every real number c L > and every measure V{ on V m without atom at lj 
and that satisfies ([8]) ; the measure on 



is exchangeable and fulfills and ([7]). 

As the erosion measure e« has infinite total mass, we see that Cj must be zero whenever /ij 
has a finite total mass. In this situation, Theorem [2] is an immediate consequence of Theorem 



3 The structure of multitype fragmentations 

The purpose of this section is to describe the structure of multitype fragmentations. In the 
monotype case, dynamics of a homogeneous fragmentation are entirely determined by an erosion 
rate c > 0, which accounts for the smooth evolution of the process, and a dislocation measure v 
on the space V m of mass-partitions, which, as its name suggests, characterizes the statistics of 
the sudden dislocations. See Sections 3.1 and 3.2 in [6]. A similar description remains valid in 
the multitype situation, more precisely dynamics are then determined by a family (cj)j 6 .n 
of erosion rates and a family (vi)ie{i,...,k} °f dislocation measures on V m , where the index i refers 
of course to the type of the fragment that is eroded or dislocated. This will be achieved first 
in the setting of partitions with types of N, and then shifted to the more intuitive framework 
of mass-partitions. 

3.1 Basic definitions 

We first introduce the natural notion of homogeneous fragmentation for mass-partitions with 
types, which bears strong similarities with that of multitype branching process. Specifically, 
let X = (X(t),t > 0) be a Markov process with values in V m and cadlag sample paths. For 
every % G {1, . . . , k}, we write Pj for its distribution starting from X(0) = 1«, i.e. at the initial 
time, there is a single unit mass with type i. 

For every mass-partition with types x = (x, i) and every real number r > 0, it will be 
convenient to write 



We then introduce a sequence of independent processes Y^ 2 \ . . . such that for every n6N, 
Y^ n > is distributed as x n X under P in . For every t > 0, we write Y(t) for the rearrangement in 
the non-increasing lexicographic order of the terms of the random mass-partitions with types 
F (1) (t),F (2) (t), . . ., and denote by P* the distribution of the process Y = {Y(t),t > 0). In 
particular P^ = Pj. 

Definition 1 The process X is called a homogeneous multitype mass-fragmentation if, 

in the sense of the Markov property, the law of X started from an arbitrary state x e V m is P s . 




m 



rx = (rx, i) = ((rx x , h), (rx 2 , i 2 ), . . .) . 
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The preceding section incites us to translate Definition [T] in the setting of partitions with 
types of N. In this direction, the notion of fragmentation operator (see Definition 3.1 in j6] on 
its page 114) has a natural extension in the multitype setting. 

Specifically, consider n G Vb a partition of some block B and 7f(0 = (n( n \n G N) a sequence 
of partitions with types. We then write Frag(7r, 7P ')) for the partition with types which is 
obtained from the collection of blocks with types of the sequence of the restrictions tt^ of 

if^ to the n-th block ir n of ix for n G N, by rearrangement in the non-increasing lexicographic 
order. In other words, each block 7r n of ix is split using Tt^ n \ Note that if the block ix n is either a 
singleton or empty, then the partition with types fc^ does not depend on tt^; more precisely, 
it is always given by ((vr n , 0), (0, 0), (0, 0), . . .). 

Let II = (Yl(t),t > 0) be a Markov process with values in with cadlag sample paths. By 
a slight abuse of notation, for every % G {1, . . . , k}, we write Pj for its distribution starting from 

n(o) = i N)i . 

Definition 2 The process II is called a homogeneous multitype fragmentation if for every 
time t > and every type i G {1, . . . , k} the distribution ofIL(t) under Pj is exchangeable, and 
the semigroup of IT can be described as follows : 

Fix t, t' > and consider a partition with types n = (ir, i) where i = (ii, i%, ■ ■ .) is a sequence 
in {0, 1, . . . , k}. Let 7fO = (ix^, . . .) denote a sequence of independent exchangeable random 
partitions with types, such that for every n G N with i n ^ 0, is distributed as fi(t') under 
Pj n . When i n = ; the block ir n is either empty or a singleton; the role ofn^ has no importance 
and its law can be chosen arbitrarily. The conditional distribution ofIl(t+t') given TL(t) = (tt, i) 
is then the law o/Prag(7r, 7P')). 

Let us now explain the connexion between these two definitions. When fl is a homogeneous 
multitype fragmentation, we know from Kingman's Theorem [T] that for every t > 0, the ex- 
changeable random partition with types Tl(t) possesses asymptotic frequencies a.s. If we write 
X(t) = |n(£)p for the random multitype mass-partition obtained by reordering these asymp- 
totic frequencies in the non-increasing lexicographic order, it can be proved that the process 
X = (X(t),t > 0) is then a homogeneous multitype mass-fragmentation. Technically, the 
main difficulty is to establish that the paths of t — > |n(t)|^ are cadlag; in this direction we 
stress that this could fail if we had equipped P m with a stronger distance such as that given 
by ([2]). In the converse direction, one can rephrase the argument of Berestycki [3] and show 
that given a homogeneous multitype mass- fragmentation X = (X(t),t > 0), there exists a 
homogeneous multitype fragmentation II such that the process (\U(t)\ l ,t > 0) is distributed 
as X. In short, there a bijective correspondence between the laws of homogeneous multitype 
mass-fragmentations and laws of homogeneous multitype fragmentations. 

Of course, the fundamental difference with the monotype case (see for instance Definition 
3.2 in [6] on its page 119) is that the distribution of the sequence tt^ which is used to split the 
partition with types n(i) into finer the blocks depends on Tl(t). However, this dependence only 
arises through the types of the blocks of U-(t), and does not involve directly the partition H(t). 
This preserves the possibility of adapting the approach developed in Chapter 3 of [6], provided 
that one can handle some technical issues. 
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In particular, it is easily seen that the fragmentation operation is compatible with the re- 
striction of partitions with types, in the sense that for every integer n 

Frag(7r,7r ( - ) )| W = Prag(7T| M , #>) . (10) 

This entails that the Markov property still holds for the restricted process ft|[ n ] = (Tl\[ n ](t),t > 
0), and as the latter can only take finitely many values, II|[ n ] is a Markov chain in continuous 
times. Note that 5|[ n ] coincides with the restriction of S|[ n+ i] to [n], and that the initial process 
IT can be recovered from the sequence of Markov chains Sim, neN. 

Just as in Section 3.1.2 of [6], these observations enable us to characterize the law of IT by 
a finite family of measures (fM)ie{i,...,k} on Vm as follows. For every fixed i £ {1, . . . , k}, n > 2 
and every partition with types 7 of [n] with 7 7^ l[n],i> we introduce the jump rate of IT|[ n ] from 

l[n],t to 7, 

By the very same arguments as in Section 3.1.2 of [6j, one can check that the collection of 
those jump rates characterize the evolution of the restricted Markov chains II|[ n i, and thus of 
the process IT. Further these jump rates can be represented as 

<?n,7 = ({* £ Vn ■ 7f|[n] = 7}) , (11) 

where /ij is an exchangeable measure on with 

jUi({l N ,i}) = and /ii ({ti £ Vn : 7f|[ n ] ^ l[n],i}) < oo for any n > 2 , 

and these requirements determine the measure fii uniquely. Note that when the condition above 
is fulfilled for n = 2, then, thanks to the exchangeability, it is fulfilled for every n > 2, therefore 
it is equivalent to (jSJ). We shall refer to the family (/ii)ie{i,...,fc} as the splitting rates of n. 

3.2 Poissonian constructions 

Our first goal in this section is to show that any family of exchangeable measures (fii)ie{i,...,k] 
which fulfill ([6]) can be viewed as the splitting rates of a homogeneous multitype fragmentation 
IT. More precisely, we shall briefly present a Poissonian construction of n which mimics that in 
Section 3.1.3 of {Bj in the monotype case. For the sake of simplicity, we assume that the initial 
state has been chosen equal to l^,^ for some type i$ £ {1, . . . , k}. 

For every type i £ {1, . . . , k}, consider the atoms (t i>m , ^ %,m \H-i y m)mm of a Poisson random 
measure in R + x x N with intensity dt <g> /Xj <g> where # stands for the counting measure on 
N. This means that for every measurable set A C R + x x N, the cardinal of the collection 
of indices m for which (ti im , Ti^' m \ £i <m ) £ A has the Poisson distribution with parameter dt <g> 
Hi®#(A) and to disjoint sets correspond independent Poisson variables. We assume that these 
Poisson random measures are independent for different types i £ {1, . . . , k}. 

For every integer n, we can construct a Markov chain in continuous time U\ n ^ = (tl^(t),t > 
0) with values in V[ n ] as follows. The atoms (t, n , £) of the Poisson point measures such that 
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7f|[ n ] = l[ n ],i or £ > n play no role in the construction and can thus be removed. Thanks to (jSJ), 
the instants t at which an atom (t, tt, £) that has not been removed arises, form a discrete set 
of M + , and the chain LT^ can only jump at such times. More precisely, if (ti tm ,TX^ l,m \i^ rn ) is 
an atom which has not been removed, then we look at the £i >m -th block of W- n ^(ti >m — ), say B. 
If the type of this block is different from the type i of the atom (in particular if B is either 
empty or a singleton), then we decide that U\ n ^(ti >m ) = U.^(ti >m —) and t ijm is not a jump time 
for the chain U\ n K If the type of B is the same as the type i of the atom, then U\ n '(ti im ) is the 
partition obtained from tv- n '(ti im —) by replacing B, that is the £i, m -th block of f[^(ti tm —), by 
the restriction of 7f(*' m ) to this block, and leaving the other blocks and types unchanged. 

To give an example, take for instance n = 7, £^ m = 2, 

vr^) = (({l,3,5,7,...},3),({2,4,6,...},l),...) , 

and set for simplicity t = t^ m . Assume also that 

= (({1, 2}, 1), ({3, 4, 5}, i), ({6, 7}, 2), (0, 0), . . .) . 

As £ i;m = 2, we look at the 2nd block of ft^ (t—), which is B = {3, 4, 5} and has type i, and thus 
coincides with the type the atom (tj, m , 7f^' m - ) , £i, m ). At time t, we split B using the partition 
with types n^' m \ This produces two new blocks with types: {3,5} which has type 3, and {4} 
which is a singleton and thus has type 0. We conclude that 

fiW(t) = (({1, 2}, 1), ({3, 5}, 3), ({4}, 0), ({6, 7}, 2), (0, 0), . . .) . 

It is easily seen that this construction is compatible with the restriction, in the sense that 
for every n G N, 5^ coincides with the restriction of I[t n+1 ] to [n]. We refer to Lemma 3.3 in 
[B] on its page 118 for the argument in the monotype case. This implies that there exists a 
process ft with values in such that the restriction of n to [n] is fl^; see Lemma 2.5 in [B] 
on its page 95 for a closely related argument. 

A crucial step is to show that for every t > 0, the distribution of n(t) is exchangeable. The 
proof relies on the following technical lemma, which is the multitype version of Lemma 3.2 in 
[6] on its page 116. 

Lemma 1 Let tx = (tt, i) G Vn be an exchangeable random partition with types and tt^ = 
G N) a sequence of random partitions with types. Suppose that : 

• it and are independent conditionally on i, 

• the sequence tt^ is doubly-exchangeable, in the sense that for every finite permutation a 
ofN, the sequences 

((T(7r (n) ),n G N) and (^ ff(n)) ,neN) 

both have the same law as if^\ Then the random partitions with types 7f and Frag(7r, n^) are 
jointly exchangeable, that is their joint distribution is invariant by the action of permutations. 

Sketch of the proof: One observes that with probability one, the conditional distribution 
of 7r given i is an exchangeable probability measure on Vfq. One can then follow the argument 
of the proof of Lemma 3.2 in □ 
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It is then easy to verify from standard properties of Poisson random measures that the 
process fl^ which has just been constructed is a Markov chain in continuous time, and that 
its jumps rates 

= fca ±P, (flW(t)=7) 

for every 7 G V[ n ] with 7 7^ 1 [■«],«, are given by ffTTj) . This shows that the process ft, which is 
specified by the requirement that its restriction to [n] coincides with W- n \ is a homogeneous mul- 
titype fragmentation with splitting rates (fH)ie{i,...,k}- Applying Theorem [2], we can summarize 
this analysis in the following statement. 

Proposition 1 Let II be a homogeneous multitype fragmentation. There exists a unique family 
( c i)ie{i,...,fc} °f nonnegative real numbers and a unique family (wi)ie{i,...,fe} of measures on V m 
which fulfill (jSj), such that the family (/ii)ie{i,...,fc} of splitting rates ofll is given by (J9j). 

Conversely, for every family (ci)j 6 /x,...,fc} °f nonnegative real numbers and every family 
( v i)ie{i,...,k} of measures on V m which fulfill (151) , if we define measures on fry ([9]), i/ten t/ie 
Poissonian construction above produces a homogeneous multitype fragmentation IT with splitting 
rates (AO»e{i,...,k}- 

Berestycki [3] established a related Poissonian construction for monotype mass- fragmentations. 
The latter can be extended to the multitype setting provided that the erosion coefficients c, are 
all the same, which enlighten the probabilistic interpretation of the dislocation measures 2^. 

Specifically, for each type i G {1, . . . , k}, consider the atoms {t^ m , x.( l ' m \ £j )m ) me N of a Poisson 
random measure in M + x V m x N with intensity di (8> z/j ® where # stands for the counting 
measure on N. Assume as usual that these Poisson measures are independent for different 
types. One can construct a pure jump process \Y(t),t > 0) in P m which jumps only at times 
ti >m at which some atom (t» )m , 5c- l > m \ £ iiTn ) occurs. The jump (i.e. the dislocation) induced by 
such an atom can be described as follows. 

We consider the mass-partition with types immediately before time t i>m , that is F(t im — ), 
and look at its £i jm -th term, say (y,j) for some y > and j G {0,...,k} (recall that the 
terms of a mass-partition with types are ranked in the non- increasing lexicographic order). If 
the type j is different from the type i of the atom, then we simply set Y(ti >m —) = Y(ti i7n ). 
Otherwise, the ^, m -th term of Y{t^ m —) is dislocated according to x.( % > m \ that is it is replaced 
by the mass-partition with types yyS t,m \ The other terms of F(tj, m — ) are left unchanged, and 
Y(ti,m) then results from the rearrangement in the non-increasing lexicographic order of all the 
terms. 

For instance, if 

4, m = 2, xfc™> = ((H,2),(i i),(o,o),...) 

and 

Y(U, m -) = f(i,4),(i i),(o,o),...^ 
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then at time tj jTn the second term of Y{t^ m — ), i.e. (|, i) is dislocated using x^ ,m ). This produces 
the sequence ((|, 2), (|, 1), (0, 0), . . .), and finally 

Y(hm)= ((\, 4),(jj,2),(i, 1),(0,0),... 

The process K is then a homogeneous multitype fragmentation with zero erosion and dis- 
location measures {vi)ie{i,...,k}- Following an argument in Berestycki [3], one can check that 
for every c > 0, the exponentially discounted process X(t) = e~ cu Y(t), t > 0, is then a ho- 
mogeneous multitype mass-fragmentation with dislocation measures iyi)ie{i,...,k} an d erosion 
coefficients q = c for every i e {1, . . . , k}. Unfortunately, this simple transformation cannot 
be extended to the case when the erosion coefficients are distinct. Informally, when the erosion 
coefficients depend on the type of the fragments, one would need information about the types 
of the ancestors of each fragment of Y(t) in order to determine the proportion of its mass that 
has been turned to dust at time t. This information is available for processes with values in 
Vfq, but not for those with values in P m . 




4 The tagged fragment 

Up to a few technical issues, the analysis of multitype fragmentations was so far an easy 
translation of that in the monotype situation. However more significant differences appear 
when dealing with finer aspects of these processes. Here, we shall focus on the evolution of the 
tagged fragment, that is the fragment which contains a point which has been picked uniformly 
at random and independently of the fragmentation process. The relevance of this study stems 
from the fact that, even though the tagged fragment alone does not characterize the evolution 
of the fragmentation, it captures some useful information. In particular, this will enable us 
to determine the asymptotic behavior of the fragmentation, by making explicit the connexion 
with multitype branching random walks. 

Let X = (X, T) be a homogeneous multitype mass-fragmentation, where X n (t) stands for 
the mass of the n- th largest fragment at time t and T n (t) for its type. It will be convenient to 
think of X as associated to a homogeneous multitype fragmentation II by X = \Tl\^. In order 
to avoid technical discussions, we shall assume throughout this section that the fragmentation 
process is conservative, i.e. for every type i G {1, . . . , k}, the erosion coefficient is Cj = and 
the dislocation measure satisfies 

Vi ^|xeP m :x = l-^x n >ojj = °- ( 12 ) 

The description of the evolution of the tagged fragment relies on the notion of Markov 
additive processes. We first provide some background in this area, referring to Section XI.2 of 
Asmussen pQ for details. 
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4.1 Background on Markov additive processes 

The class of Markov additive processes that will be useful in this work is that formed by bivariate 
Markov processes ( J t , S t )t>o, where {Jt)t>o is a continuous time Markov chain with values in 
the finite space of types {1, . . . , k}, and, roughly speaking, on every time-interval on which J 
stays constant, 5* evolves as a subordinator (i.e. an increasing process with independent and 
stationary increments) with characteristics specified by the value of J. More precisely, one 
requires that for every t, t' > 0, 

mf(St+t>) - f(S t ))g(J t+t >) | J u S t ) = E Ju0 (f(St>)g(Jt>)) (13) 

where E^- s refers to the mathematical expectation when the process (J, S) starts from the state 
(j, s) and /, g denote two generic measurable nonnegative functions. 

The law of the Markov chain (Jt)t>o is specified by its intensity matrix A = (Ajj-)jj £ {i h, 
i.e. for i ^ j, Ay is the jump rate of J from i to j and Ylj=i ^ij = 0- On everv time-interval 
[t,t + t'[ on which J = i, S evolves as a subordinator with Bernstein exponent ^W, i.e. 

E(exp(-^ i (0 )) = exp(-^ (i) (#)) . 

The Bernstein exponent is a concave increasing function which can take the value — oo, and is 
nonnegative and finite on [0, oof. 

Further, a jump of J from i to j ^ i has a probability pij of inducing a jump of 5* at the 
same time, the distribution of which is denoted by B^, and we write Bij{6) = J e^ 9x B, t j{dx). 
It is also convenient to agree that pa = 0. 

For every types i, j G {1, . . . , k} and every 9 > and t > 0, there is the following identity 
between k x k matrices : 

Me"** J* = J) = , 

where 

&(9) = -A + W) diag + (A,,/M1 " > (14) 

see Proposition 2.2 in pQ on its page 311. We shall refer to $ as the Bernstein matrix of (J, S). 

4.2 Distribution of the tagged fragment 

We are interested in the process of the asymptotic frequency and the type of the first block 
(|n 1 (t)|,t > 0) of a homogeneous multitype fragmentation n. Recall from the paintbox con- 
struction that | fl! (t) | can be viewed as the mass and type of the fragment which contains some 
point that has been picked at random according to the mass-distribution and independently of 
the fragmentation. 

The conditions which have been enforced at the beginning of this section ensure that for every 
t > 0, the first block Hi(t) is neither empty nor a singleton, hence its asymptotic frequency is 

2 Note that, since we are dealing with subordinators, we shall work with the Bernstein exponent whereas 
Asmussen pQ uses the cumulant 8 — > — 6). 
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strictly positive and its type is not 0. This allows us to introduce the process (J, S) with values 
in {l,...,fc} x R + by 

\IL 1 (t)\ = (exp(-S t ),J t ), t>0. 

Theorem 3 Suppose that the homogeneous multitype fragmentation IT has erosion coefficients 
Cj = and that its dislocation measures fulfill ({12]) . Then (J, S) is a Markov additive process 
with Bernstein matrix given for every 9 > by 



#(0) = (J_ U {i=j} - ^xl +0 l {in=j} j ^(dx 



i,je{i,...,k} 



Proof: The fact that (J, S) is a Markov process that satisfies fflBT) can be seen from the 
Poissonian construction and the arguments in Section 3.2.2 of [6]. 

The determination of the Bernstein matrix also relies on the Poissonian construction. First, 
note that for i ^ j, the jump rate Ay of the type process J coincides with the rate of occurrence 
of atoms (tj jm , 7f(*' m ), 1) with 7f( J > m ) = [n^' m \ i) and i\ = j. Using (J3J) and Theorem^ this yields 



As X/i=i = 0) this entails 



/ x nMi n =j} i'i(dx) , for i ^ j. 

J -P m \n=l / 

J V™ \n=l J 



JPm \n=] 

Thus, by (fT2"]) . we obtain the general formula 

\j = \Y1 x n*-{i»=j} ~ Mi=j} I ^( d?c ) ' (15) 

\n=l / 

A slight refinement of this argument enables us to compute the finite measure XijPijBij. 
Specifically, one finds for i ^ j 

\jPij j f{b)Bij{db) = j ^Y^x n l {in=j} f(-\nx n )^j z^(dx) . 

This gives 

X ijPij (l - ByiO)) = j (j2l {in=j} (x n - x 1 ^ 9 )) vi(dx) . 



K n=l 



Finally, the calculation of the Bernstein functions of the subordinators is made by 
reduction to the monotype situation. Specifically, we shall work under the law Pj, and we 
denote by v\ the image of Vi by the map f i : V m — > V m where fj(x) is the mass-partition given 
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by rearrangement of the terms l{j n= j}X n . Informally, this means that all the components of x 
which are not of type i are reduced to dust. Then v\ is a (monotype) dislocation measure. It 
should be plain from the Poissonian construction that if we denote by ( the instant (i.e. the 
first coordinate) of the first atom (£j jTn , Tt^' m \ 1) with 7r( l ' m ) = (ir( l ' m \ i) and %\ 7^ i, then ( is the 
first jump time of the type process J and the process killed at time (, (\Hx(t)\,u < £), can be 
viewed as the process of the tagged fragment in a homogeneous monotype fragmentation with 
no erosion and dislocation measure v\. This yields 

where, according to Theorem 3.2 in [6] on its page 135, 




^(dx) = J_ (1 - jr H^n +e ) *i(d*) • 



= /_ f E - **) i/,(dx) . 

■ / ^» \n=l / 

Putting the pieces together in (JHJ), this establishes our claim. □ 
If we introduce 



{«n=i} 



n=l 



fj(dx) < 00 for every i, j 6 {1, . . . , /c} 



then the Bernstein matrix function 6 —>■ *&(9) possesses an analytic extension to ]9_, 00 [ and 
Theorem [3] still holds for 9 e]6_, oof. 

4.3 Connexion with multitype branching random walks 

The preceding analysis of the evolution of the tagged fragment provides us with the key to shift 
some deep results on multitype branching random walks to homogeneous fragmentations. The 
approach is quite similar to that in [7], so again we shall skip details. 

Just as in the monotype case, we consider the logarithms of the masses of the fragments and 
introduce for every t > the empirical measure Z*^ = {Z^\ . . . , zj£ ), where 

00 

Z f = EVn(t)=i} 5 -ln*»(*)- ( 16 ) 
n=l 

For every fixed step-parameter a > 0, the process in discrete time (Z^ an \ n 6 Z + ) is then a 
multitype branching random walk; see [9] for a precise definition. For the sake of simplicity, we 
shall focus on the case when a = 1 and compute first a quantity of fundamental importance in 
terms of the characteristics of the fragmentation. 
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The analysis of multitype branching random walks relies on the Laplace transform of the 
intensity 

m ij (9)=E i ^e- e *zf\dx)^ = Ej ffj l {Tn(1)=j} X e n (l) 

Recall now that X(t) = \TL(t)\^ and, from that conditionally on X(t), the tagged fragment 
|5i(i)| = (exp —S t , Jt) is distributed as a size-biased sample of X(t). Hence, for every 9 > + 1, 
we have 

m l3 {0) = E, o (exp(-(0 - 1)50, = j) = (e~* ( *~% • 

We shall now make a further assumption on the fragmentation X, which will be crucial 
to investigate its asymptotic behavior. Specifically, we assume henceforth that the process J 
of the type of the tagged fragment is ergodic, i.e. the intensity matrix A = (Ay) given by 
( FToT) is irreducible. We recall from the Perron-Frobenius theory (see for instance Seneta [15] 
or Section 1.6 and II. 4 in Asmussen PQ) that for every > 0, the matrix exp(— $(0)) has 
a unique real eigenvalue with maximal modulus which can be expressed as e~ v ^ e \ In other 
words, <f (9) is the eigenvalue of the Bernstein matrix <fr with minimal real part. We also write 
u(0) = (ux(9), . . . , Uk(9)) and v(0) = {vi(9), . . . , Vk(9)) for the left and right eigenvectors El 
associated with e~^ e \ normalized so that X]i=i M «(^) = 1 an d Eti 11 '^)^^) = 1- 

We are now able to turn our attention to a fundamental family of martingales, which have 
been introduced first by Biggins in the monotype situation. 




Theorem 4 Assume that the erosion coefficients Cj are all zero, that the dislocation measures z/j 
are conservative (in the sense that f )12p holds), and that the intensity matrix (j!5j) is irreducible. 

(i) The equation 

ip(9) = (9 + l)ip'(9) 

possesses a unique solution 9 > 0. The function 9 ip{9)/{9 + 1) increases on ]9_,9[ and 
decreases on }9,oo[, and thus reaches its unique maximum at 9. 

(ii) For every 9 E]6_,6[, the process 

oo 

M(9,t)=e^Y, v ^m X n +1 (t), t>0 

n=l 

is a martingale which converges a.s. and in L l (Fi) for every type % 6 {1, . . . , k}. Further, this 
convergence is uniform for 9 in any compact set in }9_, 9[, almost surely. 



Proof: (i) It can be shown that the function 

ip :}9, oo[— > M. is concave, increasing, and ip{9) = o{9) as 9 — > oo. (17) 

It may be useful to compare our notation with that in [9j, see in particular Theorem 1 there. The matrix 
M(ff), the eigenvalue p(9) and the eigenvectors u(9) and v(9) there coincide respectively with exp(— <f>(# — 1)), 
exp(-(p(6 - 1)), u(0 - 1) and v(0 - 1) here. 
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See Theorem 3.7 in Seneta [15] for the concavity assertion, the fact that ip increases is similar. 
Finally observe from Theorem [3] that lim^oo 9~ 1 &(9) = 0, which entails <p(8) = o(9). We can 
then follow the arguments of the proof of Lemma 1 in [1] . 

(ii) We start by recall from the size-biased sampling formula (TjJ that, if (J 7 t)t>o denotes the 
natural filtration of X, then 

oo 

M(9,t) = e^X>«(*)(0)*n +1 (O 

n=l 

= e^E(exp(-eSt)vj t (6)\^ t ) . 



As v(6 l ) is a right eigenvector of exp(— 1§>(6)) corresponding to the eigenvalue e~ tv> ^ 9 \ we easily 
see from the Markov property of (J, S) that the process e tip ^ exp(—9S t )vj t (9) is a martingale 
in its own filtration. By projection on (J r t)t>o, we conclude that M(9,t) is an (jF t )-martingale. 

By an argument of discretization analogous to that in [7] , it suffices to establish the statement 
when t goes to infinity along the sequence an for some arbitrary a > 0. For the sake of simplicity, 
we shall focus on the case a = 1 without loss of generality and aim at applying Theorems 2 
and 3 of [9] to the discrete time martingale 

M(0, n) = e n ^ s) J2 v A e ) [ e- {0+1)x Z^ n \dx) , n e Z + . 

3=1 J 

Recall from Theorem 1 (ii) in [9] that Vj(9) ^ for every 9 > 9 and j = l,...,k. An 
application of the conditional Jensen's inequality to the identity 

M(6, 1) = e^ ) E(exp(-^5i)i;j 1 (^) | ^) 

shows that for every 9 > 9, there is a > 1 such that E i (M(6 l , 1) Q ) < oo for all types % — 1, . . . , k. 
On the other hand, we deduce from (i) that whenever 9 &]9, 9[, we can find a > 1 close to 1 
such that (f(9)/(9 + 1) < (f(9')/(9' + 1) where 9' = a (9 + 1) - 1. This implies that 

exp(- V (a(9 + 1) - 1) + a<p{9)) < 1 , 

and Theorems 2 and 3 in [9] now entails our claim. □ 

Let us mention an interesting consequence of Theorem H] to the rate of decay of the largest 
fragment as time goes to infinity. It follows readily from Theorem H] that a.s. 

lim -lnXi(t) = -(p'(9), 

t— >oo t 

see e.g. Corollary 1 in [5J. More precisely, the latter easily entails that a.s. for every type 
j e {l,...,k}, 

t— >oo t 

where £,•(£) = max{X n (t), T n {t) = j}. 
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4.4 Asymptotic behavior of the empirical measure 



We shall now conclude this section by presenting a couple of applications to the asymptotic 
behavior of homogeneous multitype mass-fragmentations. The first belongs to the same vein 
as Corollary 3.3 in [6] on its page 158. In the monotype case, a version of the result in discrete 
time goes back to Kolmogorov [T3], in probably the first rigorous work ever on fragmentation 
processes. Roughly speaking, Kolmogorov provided an explanation to the fact which has been 
observed experimentally in mineralogy, that the logarithms of the masses of mineral grains are 
often normally distributed. We shall see that a similar feature holds for the more general model 
of multitype fragmentations. 

Corollary 1 Assume that the erosion coefficients Cj are all zero, that the dislocation mea- 
sures Vi are conservative (in the sense that (Tl2|) holds), and that the intensity matrix (fT5|) is 
irreducible. Denote for simplicity by u(0) = u = (u\, . . . , Uk) the stationary distribution on 
{1, . . . , k} of the Markov chain J, i.e. u is the unique probability vector with 

uA = 0. 

Suppose further that if is twice differ entiable at 0. Then the following limits hold in L 2 (Pj) for 
any initial type % G {1, . . . , k} and every continuous bounded function f : R x {1, . . . , k} — > R: 

oo k 

limVl^tj/trMnl^tl.T^t)) = yW(V(0),j) 
n=l j=l 

and 

oo k 

limY,Xn(t)f(t- 1/2 (\nX n (t) + p'(Q)t),T n (t)) = Y>,E(/(jV(0, -</(0)),j)) , 

n=l j=l 

where A/"(0, — <f"{0)) denotes a centered Gaussian variable with variance —ip"(0). 



Informally, the first limit means that the masses of most fragments decay exponentially fast 
with rate f'{0) and their types are distributed according to the stationary law u of the Markov 
chain J. The second limit is a refinement of the first and shows that, in a pathwise sense, 
fluctuations are Gaussian and independent of the type. 

Proof: The two limits can be established by first and second moments estimates which rely 
respectively on the law of large numbers and the central limit theorem for the Markov additive 
process (S, J), and an argument of propagation of chaos. For the sake of conciseness, we shall 
focus on the second limit. 

Under the present assumptions, we know from Corollary 2.8 in [1] on its page 313 that as 
t — ► oo, 

*^W*=,Af(0,-^(0)), 
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where =>■ is used as a symbol for convergence in distribution. On the other hand, we also have 
J(t) =>• r, where r is a random type distributed according to the stationary law u. Further 
an easy argument using the fact that the Markov chain J mixes exponentially fast shows the 
asymptotic independence, in the sense that 

where in the right-hand side, the variables jV(0, — <f"(0) and r are assumed independent. 

Recall now that X(t) = |n(£)| J ' and, from (TjJ, that conditionally on X(t), the tagged frag- 
ment = (exp —St, Jt) is distributed as a size-biased sample of X(t). Hence 

E, ^X n (t)f(t-y 2 (lnX n (t) + <f/(0)t),T n (t)) \ = E l (j ( ^±±0)1 ^^ > 

and therefore 

(oo 
71=1 / 

= E(/(AT(0,-^(0)),r)) 

fc 

= ^ M ,E(/(AT(0,V'(0)),J))- 
j'=i 

By an argument of propagation of chaos similar to that in the proof of Corollary 3.3 in [6] 
on its page 159-160, we can estimate the second moment and get 

ME, f (jrX^fir^ilnX^ + ip'Wt),^)) 



,n=l 



X)«W(jV(0,-^(0)),j)) 

This entails the convergence in L 2 (P.;) which has been stated. □ 

Finally, using time discretization techniques similar to those in [7] , we can translate Theorem 
7 in [9] to multitype fragmentations. This yields a pathwise large deviation limit theorem for 
the empirical distribution of the fragmentation which refines considerably Corollary [lj In this 
direction, we shall further assume that the eigenvalue function ip is strictly concave and that 
the branching random walk is strongly non-lattice (see [9] for the terminology), which 
are very natural and mild conditions. Recall Theorem H] and denote the terminal value of the 
martingale M(6,t) by Af(0,oo). 



20 



Corollary 2 Let h : R x {1, . . . , k} — > R be a continuous function with compact support. Under 
the preceding assumptions, we have 

oo 

lim y/ie-W+Wm-vm y h ( t(f /(Q) + \nXJt), T n {t)) 

71=1 



V 271 ^ wi j=l 
uniformly for 6 in compact subsets of]9,9[, almost surely. 

In particular, this implies that for every a < b G R, 9 e}9, 9[ and j G {1, . . . , k}, there is the 
estimate as t — > oo 

# {n G N : ae"^' (e) < X n (t) < and T n (t) = j} 

~ ^(fl)t- 1 / 2 exp(t((0 + 1)^(0) - <p(9))) (e- a ^ - e- b ^) , 

where Ag is some strictly positive random variable with finite mean. See Corollary 2 in [9]. 
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